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1 Introduction 



Given a system of algebraic partial differential equations F = 0, where F = 
fi, . . . , fk, and a differential equation / = 0, one can effectively test if / is 
a differential algebraic consequence of F. In this paper we develop a method 
that leads to an effective procedure which finds an algebraic expression of some 
power of / in terms of the elements of F and their derivatives (or shows that 
such an expression does not exist). This procedure is called effective differential 
Nullstellensatz. A brute-force algorithm solving this problem consists of two 
steps: 

(1) find an upper bound h on the number of differentiations one needs to 
apply to F and 

(2) find an upper bound on the degrees of polynomial coefficients gi and a 
positive integer k 

such that f k is a combination of the elements of F together with the deriva- 
tives up to the order h and the coefficients gi. We solve the first problem in 
the paper. The second problem was addressed and solved in [2] and further 
analyzed and improved in |3,4J. A purely algebraic solution was given in [5]. 
Most of the references on the subject can be found in [BTTTfgfS] . 

More precisely, our problem is as follows. We are given a finite set F of differ- 
ential polynomials such that a differential polynomial / belongs to the radical 
differential ideal generated by F in the ring of differential polynomials. Know- 
ing only the orders and degrees of the elements of F and the order of /, we 
find a non-negative integer h such that / belongs to the radical of the algebraic 
ideal generated by F and its derivatives up to the order h. 

We give a complete solution to this problem using differential elimination. 
The problem is non-trivial: the first (unsuccessful) attempt was made by Sei- 
denberg [T], where it was conjectured that most likely such a bound would 
not be found. Here is where the main difficulty is coming from. In order to 
get the bound using a differential elimination algorithm we need to estimate 
how many differentiation steps this algorithm makes. Originally, termination 
proofs for such algorithms were based on the Ritt-Noetherianity of the ring of 
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differential polynomials, that is: every increasing chain of radical differential 
ideals terminates. And this result does not say when the sequence terminates. 
We overcome this problem in our paper. 

The article is organized as follows. We introduce basic notions of differential 
algebra in Section [2j Then we formulate the main result, Theorem [TJ in Sec- 
tion [31 In order to achieve this, we bound the length of increasing sequences 
of radical differential ideals appearing in our differential elimination in Sec- 
tion [5] (see Proposition [9]). For that, in Section @J we first bound the length 
of dicksonian sequences of tuples of natural numbers with restricted growth 
of the maximal element in these tuples (Lemma [8]). Finally, we apply this 
to obtain the bound for the differential Nullstellensatz in Theorem [151 from 
which Theorem [T] follows. We conclude by giving in Section [HI an alternative 
non-constructive proof of existence of the bound, based on model theory. 

There is some previous work on bounding orders in differential elimination 
algorithms. In the ordinary case, we can bound the orders of derivatives of 
the output and all intermediate steps of differential elimination [TU], and this 
bound holds for any ranking. Also in the ordinary case, one can give bounds 
for quantifier elimination [11] and for the orders and degrees of resolvents of 
prime differential ideals of a certain type [12] . A related bound for involu- 
tive prolongation, based on the analysis of stability of Spencer sequences, is 
obtained in [13J. 

Note that, unlike the bounds for differential elimination mentioned above, the 
bound for the differential Nullstellensatz proposed in this paper holds for the 
PDE case. Our bound is also based on the analysis of differential elimination. 
But, due to the ranking-independent nature of the differential Nullstellensatz, 
we could restrict our analysis to orderly rankings, which allowed us to treat 
not only the ordinary case, but the PDE well. 



2 Basic differential algebra 



One can find recent tutorials on the constructive theory of differential ideals 



One also refers to 



for differen- 



tial elimination theory. A differential ring is a commutative ring with unity 
endowed with a set of derivations A = {di, . . . , d m }, which commute pairwise. 
The case of A = {5} is called ordinary. Construct the multiplicative monoid 
6 = | 1 (?2 2 ' ' ' ^m m k» ^ o} °f derivative operators. Let Y — {yi, . . . , y n } be 
a set whose elements are called differential indeterminates. The elements of 



{9y | 9 & ©, y e Y} are called derivatives. Derivative oper- 

M = (WVi for all 0i,0 2 e and 



the set OY 

ators from O act on derivatives as 0i 
1 < i ^ n. 
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The ring of differential polynomials in differential indeterminates Y over a 
differential field k is a ring of commutative polynomials with coefficients in 
k in the infinite set of variables QY. This ring is denoted by k{?/i, . . . ,y n }. 
We consider the case of chark = only. An ideal I in k{y 1 , . . . ,y n } is called 
differential, if for all / E / and 5 E A, 5f E I. Let F C k{yi, . . . ,y n } be a 
set of differential polynomials. For the differential and radical differential ideal 
generated by F in k{yi, . . . , y n }, we use notations [F] and {F}, respectively. 

A ranking is a total order > on the set QY satisfying the following conditions 
for all E Q and u,v E QY: 

(f) 9u ^ u, 

(2) u^v^Ou^ 0v. 

Let u be a derivative, that is, u = 9yj for = df 1 ^ 2 • • • d^ 1 E Q and f < j ^ n. 
The order of u is defined as 

ordw = ord# = k\ + . . . + fc m . 

If / is a differential polynomial, / ^ k, then ord/ denotes the maximal order 
of derivatives appearing effectively in /. 

A ranking > is called orderly if ord u > ord v implies u > v for all derivatives 
u and v. Let a ranking > be fixed. The derivative 9yj of the highest rank ap- 
pearing in a differential polynomial / E k{y!, . . . , y n } \ k is called the leader of 
/. We denote the leader by Id / or u/. Represent / as a univariate polynomial 
in u f : 

f = i f u d f + amj- 1 + . . . + a d . 

The monomial is called the rank of / and is denoted by rk /. Extend the 
ranking relation on derivatives to ranks: -uf 1 > if either u\ > 112 or u\ = U2 
and di > di- The polynomial 1/ is called the initial of /. Apply any derivation 
5 E A to /: 

df 

Sf = -Tr-Suf + 5ifU d f + da^f 1 + ... + 8a d . 

The leader of Sf is 5uj and the initial of 5f is called the separantoi /, denoted 
S/. If0e0\{l}, then Of is called a proper derivative of /. Note that the 
initial of any proper derivative of / is equal to S/. 

We say that a differential polynomial / is partially reduced w.r.t. g if no proper 
derivative of u 9 appears in /. A differential polynomial / is algebraically re- 
duced w.r.t. g if deg U9 / < deg U9 g. A differential polynomial / is reduced w.r.t. 
a differential polynomial g if / is partially and algebraically reduced w.r.t. g. 
Consider any subset A C k{y 1; . . . ,y n } \ k. We say that A is autoreduced 
(respectively, algebraically autoreduced) if each element of A is reduced (re- 
spectively, algebraically reduced) w.r.t. all the others. 
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Every autoreduced set is finite [HI Chapter I, Section 9] (but an algebraically 
autoreduced set in a ring of differential polynomials may be infinite). For 
autoreduced sets we use capital calligraphic letters A, B, C, ... and notation 
A = Ax, . . . , A p to specify the list of the elements of A arranged in order of 
increasing rank. We denote the sets of initials and separants of elements of A 
by 1.4 and s^, respectively. Let = 1.4 U S4. For a finite set S of differential 
polynomials denote by S°° the multiplicative set containing 1 and generated 
by S. Let I be an ideal in a commutative ring R. The saturated ideal I : S°° 
is defined as {a G R | 3s G S°° : sa G /}. If J is a differential ideal then / : 5*°° 
is also a differential ideal (see [IS]). 

Let A = Ax, ■ ■ ■ ,A r and B = B x , . . . , B s be (algebraically) autoreduced sets. 
We say that A has lower rank than B if 

• there exists k ^ min(r, s) such that rk Ai = rk Bi for 1 ^ i < k, and 
rkA k < rkB k , 

• or if r > s and rk A; = rk Bi for 1 ^ i ^ s. 

We say that ikA = rkB if r = s and rk Ai = rk Bi for 1 ^ % ^ r. Let u be 
a derivative in k{yi, . . . ,y n }. Denote by A v the set of the elements of A and 
their derivatives that have a leader ranking strictly lower than v. A set A is 
called coherent if whenever A,B & A are such that and u# have a common 
derivative: v = ipUA = </>ub, then sbi/jA — sa4>B G (Av) : 



3 Main result 

For a finite set of differential polynomials F C k-Jj/x, . . . , y n } let -D(-F) be the 
maximal total degree of a polynomial in F . For each i, 1 ^ % ^ n, let 



For /i G Z^o let F^ 71 ^ denote the set of derivatives of the elements of F of 
order less than or equal to h. The Ackermann function appearing in our main 
result is defined as follows [271 Section 2.5.5]: 



Theorem 1 Let F C k{yx, ■ ■ ■ ,y n } be a finite set, 0^/6 {F} and let 



hi (F) — ord K (F) , H(F) 



max hi(F). 



A(0,n) 
A(m + 1,0) 
A(m + l,n + 1) 



n + 1 



A(m, 1) 

A(m, Aim + 1, n)). 
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t(F, f) be the minimal non-negative integer such that f G y TTien 
t(F, /) < A(m + 8, max(n, #(F U /), D(F U /))). 



PROOF. This result will be proved step-by-step in the following sections as 
described in the introduction and finally established in Theorem [15j 

Remark 2 It is our own choice here to bound t(F, f) using solely the maxi- 
mal orders and degrees of F and f. One might come up with another bound 
using more information of F and f. But we emphasise that the bound on or- 
ders must depend on the degrees, number of differential indeterminates, and 
number of basic differentiations as the following examples show. 



Example 3 Let f = 1 and F = jy' — l,y k \ in k{y} ; the ordinary case. In 
order to express 1 in terms of the elements of F , one has to differentiate y k k 
times. 

In the linear and non-linear cases, consider the following examples showing 
that the bound must depend on the number of variables and derivations. 

Example 4 Let f = 1 and F = {y[, yi - y' 2 , ■ ■ ■ , y n -i - Vn, Vn ~ a} in the 
ordinary differential ring k{yi, . . . , y n }, where a e k is such that a^ n ' = 1. 
We have to differentiate the first n — 1 generators n times to get y^ into the 
corresponding algebraic ideal. Hence, t(F, f) = n. 



Example 5 Let f = 1 and F = {yf,yi - yl,...,y n -i - y%, 1 - y' n } C 
k{y 1; . . . , y n }, again in the ordinary case. One can show that 

F c(y l ,y 2 ,...,y n A-y'n) = J o> 




F^d(h,yl...,y'Uy:_ 1 -2^)=L 
F^ 3) c(l 2 , y r,...,y^ 2 ,yt 1 J n 4) )=l3, 



(5) 



2) 




n-1 




2 



in — k — 1 




2»- 1 -i,2/i 



(2™- 1 ) 



n 



,V2 



(2™- 1 ) 



, . . . ,y n 



fc=i 



F <(2»-1) c (j. 



2 n -2i Vl 



(2 n -l) (2 n -l) 
1 ,V2 
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Therefore, 1 ^ (f^( 2 ' 1 1 )J. Thus, t(F, f) = 2 n , because modulo l — y' a we have 

(yTf l) = T ((y^) y' n f n - X) = T (y^f^ - ... - (2»)!(jwO' es T\y 

Example 6 If we replace F in the previous example by 

G = {u 2 xv u Xl - u 2 X2 , . . . , u Xm _ x - u 2 Xm7 1 - u x 2 n } c k{u} 

with partial derivatives d xi , . . . ,d Xm , we obtain an example which shows that 
the bound on orders must depend on the number m of derivations. Again, the 
generators will have to be differentiated 2 m times to express 1 . 



4 Bounds on lengths of sequences 

The results of this section with be further used in Section to bound lengths 
of decreasing sequences of autoreduced sets appearing in the differential elim- 
ination algorithm that we use. In this section the letters m and n will not 
mean the number of derivations and differential indeterminates, respectively. 

We begin by bounding the length of certain sequences of non-negative n-tuples. 
Call a sequence ti , ■ ■ ■ ,t k of n-tuples dicksonian, if for all 1 ^ i < j ^ k, 
there does not exist a non- negative n-tuple t such that U + t = tj. For example, 
any lexicographically decreasing sequence is dicksonian. By Dickson's Lemma, 
every dicksonian sequence is finite. Our goal is to obtain an explicit upper 
bound for the length of a dicksonian sequence, whose elements do not grow 
faster than a given function, in terms of this function, the first element, and 
the size n of the tuples. Let 

{a\, . . . ,al),{a\, . . . ,al), . . . ,{a\, . . . ,a k n ) 

be a dicksonian sequence of n-tuples of non-negative integers such that 

max (ai, . . . X) ^ fU) i 1 ) 

for all j, 1 ^ j ^ k, where 

is a fixed function. We say that the growth of this sequence is bounded by the 
function f. 

The following proposition closely resembles a particular case of our problem, 
namely that of f(i) = m + % — 1. However, in Proposition [7] the maximal 
coordinate must increase by 1 at each step, whereas in our case it is allowed to 
decrease or remain the same. We will reduce the case of a dicksonian sequence 
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with the growth bounded by a function / from a certain large class of functions 
that "do not grow too fast" , to the one treated in Proposition [7J 



Proposition 7 I28\ Proposition 1] Let t\,t2, ■ ■ ■ ,tk be a dicksonian sequence 
of n-tuples, such that the maximal coordinate of ti equals m + % — 1, for all 
1 ^ i ^ k. Then the maximal coordinate in the last tuple, t^, does not exceed 
A(n,m — 1) — 1, and there exists such a dicksonian sequence for which this 
bound is reached. 

Note that in Proposition [7J we have: m is the maximal coordinate of ti and 
the length k of the sequence is bounded by A(n, m — 1) — m. The general case 
(of any function /, not necessarily from our class) has also been studied in 
[29] using a different approach. It is shown that the maximal possible length 
is primitive recursive in / and recursive, but not primitive recursive (if / 
increases at least linearly), in n. Sequences yielding the maximal possible 
length are constructed. Moreover, if / is linear, an explicit expression for the 
maximal length is given in terms of a generalized Ackermann function. Our 
statement was motivated by the need to obtain an explicit expression for the 
bound for a wider class of growth functions. 

Let Lf jtl denote the maximal length of a dicksonian sequence of n-tuples, 
whose growth is bounded by /. For an increasing function / : Z^ — > Z^ , let 
\f^ 1 (x)] be the least number k such that f(k) ^ x. 

Lemma 8 Let f : Z^o — * Z^o be an increasing function, and d G Z^o be a 
number such that f(i + 1) — f(i) ^ A(d, f(i) — 1) for all i > 0. Then 



J f,n 



< 



f-^Ain + dJXl)-!))} (2) 



and the maximal entry of the last n-tuple does not exceed A(n + d, /(l) — 1). 



PROOF. Consider a disckonian sequence 

(a},..., 4), (al,...,a 2 n ),..., (af, . . . , a*), (3) 

whose growth is bounded by /. Construct from Q a new sequence satisfying 
the conditions of Proposition [7J Append to the first tuple d new coordinates, 
each equal to /(l), obtaining the following (n + <i)-tuple: 

(oJ,...X. /(I). •••./(!))■ 

Then add f(2) — /(l) — 1 new (n + d)-tuples. The first n coordinates of these 
tuples are (a\, . . . , a*). The last d coordinates form a dicksonian sequence of d- 
tuples, starting with (/(l), • • • , /(l)), with the maximum coordinate growing 
exactly by 1 at each step. From Proposition [7J and condition f{2) — f{\) ^ 
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A(d,f(l) — 1), such sequence exists. The last tuple will have the maximum 
coordinate equal to f{2) — 1. Next, add the tuple 



,</(2),...,/(2)). 



Since the growth of ([3]) is bounded by /, the maximal coordinate in this tuple 
equals f(2). Continue by adding /(3) — f(2) — 1 new {n + d)-tuples, whose 
first n coordinates are (af, . . . , a>D and last d coordinates form a dicksonian 
sequence growing by 1 at each step. Finally, when the tuple 



(at 



a k n J(k),...,f(k)) 



is reached, stop. We obtain a sequence of (n + <i)-tuples in which the maxi- 
mal coordinate grows by 1 at each step. We will show that this sequence is 
dicksonian. Suppose that it is not. Let tj, ti, j < I, be two (n + (i)-tuples from 
this sequence, for which there exists a tuple t such that ti = tj + t. Let t 1 , t n 
denote the first n coordinates and the last d coordinates of an (n + (i)-tuple 
t, respectively. Then we have tj —tj + t 1 and tj 
been added after the same tuple of the form 



n 



ty + t». 



If tn and ti have 



Pi 



c4/(i),...,/(i)), 



or if tj coincides with such a tuple pi and ti has been added after p i: the 
equality tj 1 = tj 1 + 1 11 contradicts the fact that the last d coordinates of the 
tuples between and Pi+i, including pi and excluding Pj+i, form a dickso- 
nian sequence. If tj and ti have been added after different tuples pi and pi>, 
the equality tj — tj + 1 1 contradicts the fact that sequence ([3]) is dicksonian. 
Therefore, our assumption was false and the constructed sequence is dickso- 
nian. 

By Proposition [71 the maximum coordinate of its last element does not exceed 
A(n + d,m — 1) — 1. Since the maximum coordinate in the first element is 
/(l) and grows by 1 at each step, the number of elements in the constructed 
sequence does not exceed A{n + d, /(l) — 1) — /(l). On the other hand, the 
number of elements in the constructed sequence is: 

/(2) - /(l) + /(3) - /(2) + . . . + f(k) -f(k-l) + l = f(k) - /(l) + 1. 

Therefore, 



that is, 
and 



/(*) - f(l) + l^A(n + d, /(l) - 1) - /(l), 
f(k)<A(n + dJ(l)-l), 



k < 



>- 1 (A(n + d,/(l)-l 
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5 Differential elimination algorithm 



Using the result of the previous section, we obtain an upper bound for the 
length of sequences of autoreduced sets of decreasing rank produced by a 
differential elimination algorithm. The idea is to put in correspondence with 
such a sequence a dicksonian sequence of tuples, whose growth is bounded by 
a function derived from the algorithm. 



We fix an orderly ranking. Algorithm [T] computes a characteristic decompo- 
sition of a radical differential ideal given by a set of generators. It is designed 
in such a way that allows us to control the orders and degrees of differential 
polynomials occurring in all intermediate steps together with a bound on the 
number of iterations of this algorithm. Also, in the algorithm the procedure 
algrem computes an algebraic pseudo-remainder of a polynomial with respect 
to an algebraic triangular set (a set is called triangular if the leaders of its 
elements are distinct). Algorithm MinimalTriangularSubset inputs a finite set 
of differential polynomials and outputs one of its least rank triangular subsets. 
Algorithm CharSet inputs a finite set of differential polynomials and outputs 
one of its characteristic sets, that is, an autoreduced subset of the least rank. 
Denote by A(C) the set of "differential S-polynomials" of C defined in [T4l Def- 
inition 4.2]. From now on, m and n again denote the numbers of derivations 
and differential indeterminates, respectively. 



Algorithm 1 RGBound(Fi) 
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Input: a set F x C k{y 1; . . . , y n } with derivations {di, . . . , d m } 

Output: A finite set T of triangular sets such that {Fi} = f| [C\ : ; 

cer 

if 1 ^ [C] : then C is coherent and autoreduced 

otherwise 1 G (C) : 
T:={0}; ^:={(F!,0)} 
while [/ 7^ do 

Take and remove any (F, C) G U 

f := an element of F reduced w.r.t. C of the least rank 
if s f $ k then U := U U (F U S/, C) end if 
if i 7 ^ k then Z7 :=[/U(FU i /; C) end if 
D:={C eC\ ldC= Id/ for some G 0} 
C:=C\DU{f}; G:=FUA(C)UD\{/} 
6 := maxordg 

B := MinimalTriangularSubset({^C | C G C, oidOC < ft}) 
B:= {algrem(/i,£\{/i}) | h G i3} 
if rk£ ^ rk£ then T := T U {£}; continue; end if 
F := {algrem(#, 6)| 5 gG}\ {0}; C := CharSet(s) 
if 72 = then T := T U {C} else t/ := [/ U (F U F, C) end if 
end while 
return T 

We get the following bounds for the growth of the maximal degrees of the 
polynomials computed at the z-th step of Algorithm [TJ 

Proposition 9 Fix (F;,Cj) ^ G U and let (F + i,Cj + i) ^ be any of the 
elements obtained from (Fj,Cj) after one iteration of the while-Zoop. We then 
have 

D{F i+1 UC i+1 ) <(4D(FiUCi))( ™ 

PROOF. Consider the iteration of the loop. In the first five lines of the loop 
the degrees do not change as adding an initial or a separant does not cause an 
increase in the degrees. So, the first place where the degrees may change is the 
computation of A (C^j , that is, computing cross differentiations. This at most 
doubles the degrees. After that, the only places where the degrees of polyno- 
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mials may increase are calls to algrem(g, B) for g G G or algrem (h, B \ {h}) for 
h G B. In both cases it is a sequence of at most \B\ algebraic pseudodivision. 
We will prove the bound for the reduction of a fixed g e G modulo B, and the 
other case follows similarly. 

Assume that \B\ = N, and let B = {Bi, . . . , B^} be ordered such that 

ld(-Bi) > \d(B 2 ) >■■■> \d(B N ). 

Let g^ := g and g^> := algrem (g, {B\, . . . B t }) for t > 0. Denote the maximal 
total degree of g® U B by 5(t), t ^ 0. Note that 5(0) «C 2D(Fi U C,). Then 
j s obtained by the pseudo-division of with respect to the polynomial 
B t+1 . Thus, 

9 {t+1) = ^ +1 9 {t) - QB t+1 , 

where e is a sufficiently large exponent specified below, q is the pseudo- 
quotient, and the degree of in ld(I? t+ i) is smaller that the same degree 
of £>t+i- The exponent e is bounded by 



degid( Bt+l) (g (t) ) ~ deg ld(Bt+l) {B t+1 ) + 1 s= 5{t). 
Therefore, the total degree of is bounded by 



8{t + 1) < 6(t)6(0) + 5(t) + 5(0) < 



f 25{0)5{t), t^lox 5(0) ^ 2; 
[3, t = and 5(0) = 1. 



This implies that 5(N) ^ (25(0)) N+1 . Using that 5(0) ^ 2D(F t U &), N ^ 
( fei ^ m ) ; an d h ^ 2H(Fi U Cj), we get the claim (the numbers bi are defined in 
Algorithm [T]). 



5. 1 Differential bounds for splitting 



Consider now the splitting part of differential elimination. Algorithm [T] re- 
moves an element (F, C) from U and within one iteration of the while-loop it 
converts this element into one, two, or three elements. Moreover, if the set C 
does not change, the orders and degrees of the elements of F do not increase 
after the conversion. We call such an iteration incomplete. Denote the maxi- 
mal number of all (complete and incomplete) iterations of the while-loop of 
Algorithm [1] by L(F). 

Proposition 10 Algorithm^ is correct and terminates. Moreover, 

L(F)<:iog 2 (A(m + 7,Q(F)-l)), 
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where 



Q(F) = max ^9, n, 2 9H( - F \ D(F)j . (4) 

PROOF. To demonstrate correctness we will show that the while-loop of 
Algorithm [1] has the following invariant 

{F 1 }= n {FUC}:H? f)[A}:H%. (5) 

(F,C)eU AdT 

Indeed, since F% C F for all (F, C) G U and for any A G T every element of 
F\ is reducible to zero with respect to A, we have the inclusion "c" . We will 
show the opposite inclusion by induction on the number of iterations of the 
loop (not assuming that it is finite). Invariant ([SD holds at the beginning of 
the first iteration (T — 0,U = (F\,0)). Let a finite number of iterations of 
the loop be executed preserving the invariant. Suppose that at the next step 
we remove an element (F, C) from U. Let C be either the set added to T (in 
this case we let F = 0), or the second element of the pair {F,C^j returned 
back to U. Note that in both cases we have 

{F,C}C{F,C} (6) 

by construction. We will show the inclusion 

{F, C) : Hf H {F, i f ,C}: ^ c °° n {F, s f , C} : C {F, C} : H™, 

that together with the inductive hypothesis shows the result. Indeed, apply- 
ing [HI Proposition 6.6], since the polynomial / chosen in the loop is reduced 
with respect to C, we have 

{F, C} : (H c U H f )°° fl{FU i f , C} : R {F U s f , C} : = {F, C} : . 

It remains to note that 

{F,C} : Rf C {F,C} : {H c UH f }°°. 

Indeed, according to Algorithm [TJ every g G C comes from some C G C as a 
remainder with respect to the triangular set Be [C]. Applying (jSj) together 
with (TUJ Lemma 5] and [HJ Lemma 6.9] to K = s g , H = He, we obtain that 

{F, C) : H C C {P, C) : [H s U F c )°° C 

C {F, C} : (He U ifc) 00 = {i 71 , C} : (F c U Hf)°°. 

It remains to note that in the case of rk B 7^ rk B we have 1 G (£>) : by [TUJ 
Lemma 5]. 
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Termination of the algorithm and the bound for L(F) will be proved as 
follows. To each element of U we associate an (m + 4)-tuple in such a way that 
the sequence of these vectors for each element of U is dicksonian. Note that 
by our definition L(F) is the maximal length of such dicksonian sequence. At 
the beginning, the set U consists of one element only. We associate the vector 

r = (0,...,0,n,n,D(F)) 

to it. Let an element (F, C) with vector r be under processing of the loop and 
let / be an element of F reduced with respect to C and of minimal rank with 
this property. 

If an incomplete iteration occurs, we add an element of the form (FU{g}, C) to 
U. The degree of g is less than the one of / (this is either the initial or separant 
of /) . We then transform the vector r by replacing the last its component by 
D(g). If the iteration is complete then U is concatenated with a converted 
element (^F,Cj. In this case we change the components of r in the following 

way. Let rkf = {Oyj) d be the rank of /, where 9 = c^ 1 . . . d%?. Consider the 
(m + 4)-tuple 

f = (i u ... } i m} d } j } n-j, deg(g)), 

where g is the element of F reduced with respect to C and of minimal rank. 
Proceeding in the described way from r to f, let r k denote the (m + 4)-tuple 
corresponding to the fc-th iteration of the while-loop. Since / is reduced w.r.t. 
C, the sequence r , r 1; r 2 , . . . is dicksonian. Therefore, the number of iterations 
for each element is bounded. This proves termination of Algorithm [TJ Note 
that if we remove n — j from f, the sequence might not be dicksonian. Indeed, 
let F = y 2 } C k{yi,y 2 ,y 3 } with y 1 < y 2 < y 3 . We then would have 
n = (0,1, 1,1) and r 2 = (0,1, 2,1). 

Let Hi = max(F(F 1 ),m), H k+1 = 2H k , D x = D(F 1 ), and 

By Proposition [91 the maximal coordinate of the [m + 4)-tuple r k does not 
exceed max(H k , D k ,n). Let 

Ui = max{n,9,2 9H{F \D{F)) , u k+1 = 2^ 2+lo ^ Uk \ 

Then the maximal coordinate of r k does not exceed u k for all k ^ 1. Indeed, 
we will prove by induction that H k ^ | log 2 u k and D k ^ u k . For k = 1 these 
inequalities hold by definition of u\. Assuming that they hold for H k , D k , and 
Uk, prove them for H k+ i, D k+ i, u k+ \. Since Hk+i = 2H k , we have 

2 9H fc+1 = 2 9-2H k = (^H^ 2 ^ ^2 = 2 2\o g2 u k ^ 2 ^log 2 « fc < 
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because u k ^ 9. Next, 

' ( 2H k + m> 



log 2 D k+ i 



m 



+ 1 (2 + log 2 D k ) <: 2 2H * +m (2 + log 2 u k ) ^ 



< 2 3H *(2 + log 2 u k ) ^ ^(2 + log 2 u k ) = log 2 ttfc+i. 

Here we used the fact that Fffc ^ m, as well as the inequality H k ^ | log 2 
proven above. Now observe that for x ^ 9, we have 

2 ^(2+io g2 x) ^ 2 x+2 _ 3 = ^4( 3):c _ !) 

Therefore, the sequence of (m + 4)-tuples To, t±, r 2 , . . . satisfies the conditions 
of Lemma [8] with <i = 3. And, according to this lemma, the length of this 
sequence does not exceed 

f-HA(m + 7, 



where f(k) = u k ^ 2 k . We can now plug in /(l) = u 1; and replace / 1 with 
log 2 - 

Corollary 11 The maximal orders and degrees of polynomials computed by 
Algorithm^ do not exceed 



A(m + 7,Q(F) - I) 



(7) 



PROOF. Follows directly from Lemma [Hand Proposition [TU1 



5.2 Lifting the final bound from splitting 



Note that for / G k{yi, . . . , y n } and F C k{yi, . . . , y n } we have 

1 G {F} : / / G {F}. 



Lemma 12 Lei Ai, . . . , A p = .4. C [F] 6e a coherent autoreduced set that 
reduces all elements of F to zero. Then for all f G {F} 

1 G (A™) : m U /) , 
where q := max I 0, ord/ — minord g ) . 
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PROOF. By our assumption, 1 G {F} : /. Moreover, since [F] C [A] : H% , 
it follows from [T9J Section 5.2] that 

{F}:f= [A] : (H% U /) n fU^} : / n f]{F, Si } : /, (8) 

i=l i=l 

where ij and Sj are the initial and separant of A», respectively. Hence, 1 G 
[y4] : U /), that is, / G [A] : H°g- Let g be a partial pseudo- remainder of 
/ with respect to A. Then by the Rosenfeld lemma (THJ Lemma 5, III. 8], we 
have g G (*4) : HJ£- Since the ranking is orderly, there exsits h G such 
that h- f — g G L4^ 9 M. Indeed, at each step of the partial pseudo-division the 
order of the resulting differential polynomial is less than or equal to the one 
of the previous polynomial. And q represents the maximal number of times 
one possibly needs to differentiate elements of A to perform one step of the 
partial psuedo-reduction. 

Lemma 13 Let F C k{y 1 , y n ] and f, f u . . . , f k G k{y u ...,y n }. Suppose 
that for some d G Z^i we have (fx ■ . . . ■ ft) d G (F) : f . Then 

9ih-...-e k f k e^(F(^ k+1)H+ld )):f, 
where 9i G O with ord#j ^ H for all i, 1 ^ i ^ k. 

PROOF. It follows from the proof of [301 Lemma 1.7] that if a d G (F) then 
for any a G k{?/i, . . . , y n } and 9j G A. Therefore, 

mh) ■ h ■ ■ ■ ■ ■ fk + ■ ■ ■ + fi ■ ■ ■ ■ ■ fk-i ■ (ft/*)) 2 " -1 G (F^) : f°°. 
Multiplying by ((dji) ■ f 2 ■ . . . • /fc) 2d ~ 1 , we obtain 

mh) ■ h ■ ■ ■ ■ ■ fkf {2d - 1] e (f^) : r. 

To make the computation simpler, we have 

{{d t h)-h-...-f k f d e (f™) :/°°. 
By induction, we conclude that 

mh) ■ h ■ ■ ■ ■ ■ hf Hd e (V(<< i+4+ - +4H ))) : r c (V(^ 4H+ld )) : r. 

Similarly, we obtain 

(PiA) • {0 2 f 2 ) ..... / fc ) 4 ^ G ( F («** +l ^* +1 4**)) : /oo = 
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Finally, by induction we get 

((fc/O • (6 2 f 2 ) ..... (6 k f k )f Hd e (V(^ +1 <i+4* + ... + 4(^*))^ . r c 




which finishes the proof. 

Lemma 14 For a finite subset F C k{?/i, . . . , y n } we /iai>e: 

t (F, /) < ord / + HOT ■ 2^> + 4 (-"' r >"" ,+ '«)-'< G '"« (i , (9 ) 

<* := max (£>(/), A (m + 7, Q(F) - ^n**^^-*, ; (1Q) 

G C k{yi, . . . ,y n } is such that L(G) ^ L(F) - 1 and H(G) ^ H(F) ■ 2 L ^ 
and Q(F) is defined in formula (TjJ. 

PROOF. Let A be the first component computed by Algorithm [I] and p be 
the number of elements of A. Note that 

p <: n ■ 2 H{A)+m . (11) 

We then have 

{F} = [A]:H%n(\{F,M}nf){F, Si }. 

i=i i=i 

If le(A):H%, we have 

si s p - ii i p e \[{A). 

Therefore, 

/ ■ si s p - ii i p e \f(A). 
Consider now the case when A is as in Lemma [12] that gives us 

/•si 1 -...-s>-i^-...-i^G (^ ord/) ) 

for some non-negative integers ji, . . . ,j p and fci, . . . , k p . Therefore, 

/ • Si Sp ■ u Ip e y/(A(***fi). 

Hence, by [5j Corollary 1.7], which gives an upper bound for degrees in the 
algebraic Nullstellensatz when the degrees of the generating polynomials are 
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not equal to two, after squaring all elements of „4^ ord /) of de gree two, we 
obtain that 

(/•s 1 -...-s p -i 1 ....-g d G (a^°^), (12) 

where d is defined by fflQj) . Indeed, n ■ 2 H ^ +m+ord f bounds the number of 
algebraic indeterminates in ^ ord /) and /, H(A) ^ H(F) ■ 2 L ^ (b ecause 
at each step of Algorithm [1] the maximal order doubles at most), and by 
Corollary [11] we have 

max (4, D (A^ ovAf) )) = max(4, D(A)) ^ 

^A(m + 7,Q{F)-1). 

We also have 

le{F,8i}:fmdle{F,u}:f 

for all i, 1 ^ % ^ p. Let G be F, ij or F, Sj with the maximal t(G,f). It then 
follows that 

I h = U + K 



where f 3 G y/(F<.<X G >f))), hj G J (ip (GJ)) ) or ^(sp (G,/)) ) , and kj is a natural 
number. Multiplying the above expressions, we obtain that 

f^=g + h, (13) 

where g G (f ( ^* (g,/)) ) and 

G ((0 lSl ) • . . . ■ (fl p s p ) • (fliu) • . . . • (^i p )| ord(^), ord(^) < t(G, /)) . 

Moreover, again since at each step of Algorithm [JJ the maximal order doubles 
at most, we have 

(A {siq) ) C (F^ q+H{F) - 2HF) A 
for any q G Z^o- By Lemma [TBI and inclusion (fl2l) we have 



lSl ) (fl p S p ) ■ fe) (e' p l p ) G Jf A «ord/ + 4(^)-KO,/) + X,)>| ; f 



where d is defined in (fTDl) . Hence, 

(MO ■ • • • • (0 P s p ) ■ (^n) ..... (#;g g 

! / . .1 , . ... i r. 



p(^ovdf+H(Fy2 L ( F ')+4( 2 P+ 1 )- t ( G 'f)+ 1 d)\ . y. 



18 



for all 6 k and 6' k with ord(^), ord(0£) ^ t(G,f), 1 ^ k ^ p. Thus, from 
inequality (|TT|) and inclusion (fl4l) it follows that 



h e 



\ 



^^ord/ + H(F).2MF) +4 ( 2 - 2ff( ^ )+1 )- t ( G '« +1 ( i 

V 

which finishes the proof because we have representation (1131) and g G (^F^ 1 ^^^ 

Theorem 15 We have 

t(F, f) <: A{m + 8, max(n, H(F U /), L>(F U /))). (15) 

PROOF. We begin the proof by recalling Corollary [TT], which states that at 
any stage of Algorithm [1] the orders and degrees of differential polynomials 
computed by this algorithm do not exceed the bound 

E:=A(m + 7,Q(FUf)-l), 

where the function Q is defined in formula (jl]) and Q(F U /) ^ Q(F). In 
particular, 

ord/ E, H(F) ^ E, 

and by Proposition [TU1 

L(F) ^\og 2 E. 
It follows directly from the definition of Q(F) that 

n^E, m < £, 100 < £. 

Using these inequalities, we can bound the quantity d defined in ( TTUI) as 

d < F e ' 2e2+2E ^ E eeE , 

whence 

d . 4 d ^ 4 2d < . 
To simplify the latter formula, we note that 

p times pk times 

for all natural numbers k and p, which can be easily derived by induction, 
using the inequality ah ^ a b , which holds for all natural numbers a, b ^ 2. 
Thus we get: 

d-4 d ^4 2d ^ 2^ . 

5E times 
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This allows us to obtain the following inequality from (Q: 

, \t(G,/)+l / \t(GJ)+l 

t(FJ) + l^E + l + E 2 +i 2^ \ ^ I 3 • 2^ ^ 

\5E times/ \ 5E times/ 

t(G,f)+l 



\6E times/ 

Now we use the fact that 



2 ... 



A(A,k)= -3 

fc+3 times 

and simplify the above formula as 

t(FJ) + l<:A(4,QE) t(G ' f)+1 , 
and noting that A(4, 6£) «C A(4, A(5, E-\)) = A(5, E) yields 

t(FJ) + l ^A(5,E) t( - G >n +l . 

Now recall that the set G, defined in the proof of Lemma [HI is the set obtained 
from F by Algorithm [T] by adding to it an initial or separant. If now we take 

E G :=A(m + 7,Q(G)-l), 

and let Gi be the set obtained from G at the next iteration of the Algorithm [T] 
by adding an initial or separant, we can similarly write 

*(£,/) + 1< A(5,E G ) t( > Ch >n +1 . 

We continue recursively writing similar inequalities for G 2 , G 3 , . . ., noting that 
the length of this chain does not exceed the number of iterations in Algo- 
rithmdj that is, L(F) ^ log 2 E. We note also that all quantities E, Eq, Eq 2 , ■ ■ ■ 
arising in these inequalities can be uniformly bounded by 

E := A (m + 7, max (9, n, 2 9E , e)) = A (m + 7, 2 9E ) ^ A(m + 7, A(A, E)), 

since the orders and degrees are uniformly bounded by E. Thus, we have 

t(F, f) + 1 < W^Z, 

E times 
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where W = A(5, E), which implies that 

t(F,/) + l< 2^ ^ A(A,WE) ^ A{A,W 2 ) ^ A(A,A(5,W -1)) = 

WE times 

= A(5, W) < A (5, A (6, E - l)) = A (6, = 

= A(G, A(m + 7, A(4, £))) < A(m + 6, A(m + 7, A(4, £))) = 

= A(m + 7,A(A, E) + l)) = 

= A(m + 7, A(4, A(m + 7, Q(F U /) - 1)) + 1) < 

sC A(m + 7, A(m + 6, A(m + 7, Q(F U /) - 1)) + 1) = 

= A{m + 7, A(m + 7, Q(F U /)) + 1). 

Note that 

Q(F U /) < A(4, 5) < A(m + 8, 5 - 1) - 1, 
where 5 := max(n, U /), D(F U /)) ^ 1, since n ^ 1. Therefore, 

A(m + 7,A(m + 7, Q(F U /)) + 1) < 

^ A(m + 7, A(m + 7, A(m + 8, 5 - 1) - 1) + 1) < 
< A{m + 7, A(m + 7, A(m + 8, B - 1)) - 1) = 
= A(m + 7, A(m + 8, 5) - 1) = 
= A(m + 8,5). 



6 Model-theoretic proof of existence of the bound 

The following argument was shown to the authors by Michael Singer. In this 
section we refer the reader for ultrafilters and construction of ultraproducts 
to books in model theory, for instance [3"Tf52] . Let k be the differential closure 
of k (see [551 Definition 3.2] and the references given there) and q G Z^ - We 
would like to emphasise that in the statement below we had to fix in advance 
the number r of differential polynomials in F to be able to use ultraproducts. 
So, in Theorem [16] the variable r is quantified before the bounding function 
f3. However, the constructive bound that we obtained in Theorem [T51 does not 
have such a restriction, because it depends solely on the orders and degrees in 
F and /, but not on the number of elements in F. 

Theorem 16 For every r G Zj> there exists a function (3 : Z^ 3 ^0 such 
that for any q G Z^ and F C k{y%, ■ ■ ■ , y n } with 

\F\ = r, max(H(F),D(F)) ^ q, and 1 G [F] 

we have 

1 e ( F (^l3(d,r,n))\ _ 
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PROOF. Assume that the statement is wrong, that is, there exist r, q G Z^o 
such that for any a G there exist p lia , . . . ,p r , a G k{?/i, • • • 5 2/n} with 

max(iy(p ii ),D(p^)) < g 

such that 

1 G [Pl,a, • • • ,Pr,a] 

and 

MEEwS ( 16 ) 

i=l j=0 

for all G k{?/i, . . . ,?/„,} of order less than or equal to q + a. Again, it is 
essential here that r does not depend on a. For a maximal differential ideal 
M in the differential ring IIiez> k denote the differential ring (lliez >0 kj /M 
by -R'm- There is a natural differential ring homomorphism 

II ^ ■ ■ ■ ' Vn} ^m{?/i, ■■■,!/»}=: 

We shall now make a special choice of the maximal differential ideal M. Let 
T be the filter consisting of all cofinite subsets of Z^o- Then, there exists an 
ultrafilter U containing T . Since the field k is differentially closed, by Los' 
theorem [3T1 Theorem 8.5.3] the ultraproduct 

K := Hk/U 

is a differentially closed field with the following property. 

Let a = (do, ai, ■ ■ .) and b = (bo, b\, b 2 , . . .) G K. Then, we have: if a — b 
then cij = bi for infinitely many indices i. We now take M to be the kernel of 
the differential ring homomorphism 

n k-K 

«eZ;> 

Let pj be the image of (pi ; x, Pip, Pi,s, • • •) m R- This is defined correctly as all 
Pij have orders and degrees bounded. Assume that (z%, . . . , z n ) G (K M ) n is a 
zero of pi for all i. Then, for each i, 1 ^ % ^ r, there exists VJ C Z^o, VJ G W, 
such that 

Pi,j{.^l,j , ■ ■ ■ ) Z n ,j) 

for all j G Vi, where (^,1,^,2,-^,3, • • •) is mapped to 2 f under the mentioned 
differential ring homomorphism for each t, 1 ^ t ^ n. Since Vi D . . . D V r G U 
and tf:U, there is an index 

j G Vi n . . . n V r . 
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Therefore, 

Pl,j\%l,jj • • • > 2 n j) ... "p T j{Z\j^ ■ ■ ■ i Zn,j) 0. 

Since k is differentially closed, this contradicts to 1 G \pij, ■ ■ ■ ,Pr,j] m the 
differential ring k{yi, • • • , Vn}- Therefore, since the field Km is differentially 
closed, we have f G [pi, . . . ,p r ). Hence, there exist 7 G Z^. and differential 
polynomials qij G Km{vi, ■ ■ ■ ,Un} with ord qlj < 7 + q so that 

i=l j=0 

Again, due to our choice of M (that is, due to the fact that U is an ultrafilter) , 
there exists a G Z^o with a > 7 such that 

r 7 

i=i 3=0 

where q^j G k{yi, ...,?/„} of order less than a + q. Since pj jQ , G k{yi, . . . , y n } 
for all i, 1 ^ % ^ r, by taking a basis of k over k we may assume that in fact 
qij G k{yi, . . . , y n } for all % and j, 1 ^ i ^ r, ^ j ^ 7. This contradicts 
to (fTBj) . Thus, our initial assumption was wrong. 



7 Conclusions 

We have obtained the first bound on orders for the differential Nullstellensatz. 
Surely, one can improve the bound and find many applications of it. A general 
programme which is being realized here is as follows. The differential elim- 
ination algorithms would be very useful for applications if there were faster 
versions of them. Our work on bounding orders could lead to: 

(f ) understanding complexity estimates for the differential elimination, 
(2) developing combined and separated differential and high performance 
algebraic algorithms. 

One of the ideas is, instead of using the usual differential elimination, perform 
all differentiations at the beginning of the process and then use only fast alge- 
braic methods. We hope that our bounds will contribute to this programme. 
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